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Abstract 

Quantal effects on growth of spinodal instabilities in charge asymmetric nuclear matter are 
investigated in the framework of a stochastic mean field approach. Due to quantal effects, in both 
symmetric and asymmetric matter, dominant unstable modes shift towards longer wavelengths and 
modes with wave numbers larger than the Fermi momentum are strongly suppressed. As a result of 
quantum statistical effects, in particular at lower temperatures, magnitude of density fluctuations 
grows larger than those calculated in semi-classical approximation. 
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I. INTRODUCTION 



In many processes, such as induced fission, heavy-ion fusion near barrier energies and spin- 
odal instabilities and nuclear multi-fragmentation, dynamics of density fluctuations play a 
dominant role. For description of these processes mean-field transport models, such as 
;ime-dependent Hartree-Fock (TDHF) Q] and the Boltzmann-Uhling-Uhlenbeck (BUU) 
3I models are not very useful. TDHF includes, the so called, one-body dissipation mecha- 
nism, but associated fluctuation mechanism is not incorporated into the model. Similarly, 
the extended TDHF and its semi-classical approximation BUU model involves one-body 
and collisional dissipation, but the associated fluctuation mechanisms are not included into 
the description. It is well known that no dissipation takes place without fluctuations. In 
order to describe dynamics of density fluctuations, we need to develop stochastic trans- 
port models by incorporating fluctuation mechanisms into the description. There are two 
different mechanisms for density fluctuations: (i) collisional fluctuations generated by two- 
body collisions and (ii) one-body mechanism or mean-field fluctuations. Much effort has 
been given to improve the transport description by incorporating two-body dissipation and 
fluctuation mechanisms. The resultant stochastic transport theory, known as Boltzmann- 
Langevin model 4, (|, provides a suitable framework for dynamics of density fluctuations 
in nuclear collisions around Fermi energy. However, two-body dissipation and fluctuation 
mechanisms do not play an important role at low energies. At low bombarding energies, 
mean-field fluctuations provide the dominant mechanism for fluctuations of collective nu- 
clear .notion. In a recent work, we addressed this q oesrion g Restricting on, treatment 
at low energies, we proposed a stochastic mean-field approach for nuclear dynamics, which 
incorporates one-body dissipation and fluctuation mechanisms in accordance with quan- 
tal dissipation-fluctuation theorem. Therefore, the stochastic mean-field approach provides 
a powerful microscopic tool for describing low energy nuclear processes including induced 
fission, heavy-ion fusion near barrier energies and spinodal decomposition of nuclear matter. 

Much work has been done to understand the spinodal instabilities and their connection 
with liquid gas phase transformation in symmetric and more recently charge asymmetric nu- 
clear matter. Most of these investigations have been carried out in the basis of semi-classical 
Boltzmann- Langevin (BL) type stochastic transport models [8!]. There are two major prob- 
lems with these investigations. First of all, numerical simulations of BL model are not very 
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easy, even with approximate methods, simulations require large amount of numerical effort. 
The second problem is related with the semi-classical description of spinodal decomposition 
of nuclear matter. According to our previous investigations, quantal statistical effects play 
an important role in spinodal dynamics 0, [lo|. There are qualitatively two different regimes 
during evolution of nuclear collisions in Fermi energy domain. During the initial regime 
of heavy ion-collisions, namely, from touching until formation of hot and compressed piece 
of nuclear matter, collisional dissipation and fluctuations are substantially important. On 
the other hand, during expansion of the system into mechanically unstable spinodal region, 
collisional effects may be neglected. In the spinodal region, local density fluctuations, which 
are accumulated during the initial regime, are mainly driven by the mean- field until system 
breaks up into clusters. Recently proposed stochastic mean-field approach provides a useful 
tool for describing spinodal decomposition of expanding hot piece of nuclear matter. The 
approach includes quantum statistical effects and at the same time, numerical simulations 
of the approach can be carried out without much difficulty. 

In this work, we study early growth of density fluctuations in charge asymmetric nuclear 
matter and investigate quantum statistical effects on spinodal instabilities and on growth 
rates of dominant unstable modes on the basis of stochastic mean-field approach. In sec- 
tion 2, we present a brief description of the stochastic mean-field approach. In section 3, 
we calculate early growth of density fluctuations, growth rates and phase diagram of dom- 
inant modes in charge asymmetric systems, and study quantal effects on these quantities. 
Conclusions are given in section 4. 



II. STOCHASTIC MEAN-FIELD APPROACH 



In the standard TDHF description of a many-body system, time-dependent wave function 
is assumed to be a single Slater determinant constructed with time-dependent single-particle 
wave functions. The standard approach provides a good description for the average evolution 
of collective motion, however it severely restricts fluctuations of collective motion ll, |2j. In 
order to describe fluctuations, we must give up single determinantal description and consider 
superposition of determinantal wave functions. In the stochastic mean-field description, an 
ensemble of single-particle density matrices associated with the ensemble of Slater deter- 
minants is generated in a stochastic framework by retaining only initial correlations Q. A 
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member of single-particle density matrix, indicated by label A , can be expressed as, 

In this expression and in the rest of the paper label a = n,p represents neutron and proton 
species and < i\Pa{Q)\j > are time-independent elements of density matrix determined by 
the initial correlations. The main assumption of the approach is that each matrix element 
is a Gaussian random number specified by a mean value < i\p^(0)\j > = 5ijp a (i) and a 
variance, 

< z|^(0)|j >< j'\6f${0)\V > = l -5 ab 8 iV 5 jjr { Pa {i)[l - Pa (j)} + Pa (j)[l - Pa (i)}}. (2) 

In these expressions < i\5p^(0)\j > represents fluctuating elements of initial density matrix, 
p a (j) denotes the average occupation number. At zero temperature, the average occupation 
numbers are zero and one and at finite temperature, they are given by the Fermi-Dirac 
distribution. In each event, different from the standard TDHF, time-dependent single- 
particle wave functions of neutrons and protons are determined by their own self-consistent 
mean-field according to, 

ih^(f,t;X) = h x a ^(r,t;X). (3) 

Here = p 2 /2m a + U a (n x ,rip) denotes the self-consistent mean-field Hamiltonian in the 
event, which depends on proton and neutron local densities n„(r, t). We can express stochas- 
tic mean-field evolution in terms of single-particle density matrices of neutrons and protons 

as, 

ifi§ i P X a (t) = [h x a ,p x a (t)]- (4) 
In the stochastic mean-field approach an ensemble of single-particle density matrices is 
generated associated with different events. In this approach, we can calculate, not only 
the mean value of observables, also probability distribution of observables. Even if the 
magnitude of initial fluctuations is small, in particular in the vicinity of instabilities mean- 
field evolution can enhance the fluctuations, and hence events can substantially deviate from 
one another. By projecting on a collective path, it is demonstrated that the stochastic mean- 
field approach incorporates one-body dissipation and one-body fluctuation mechanisms in 
accordance with quantal dissipation-fluctuation relation Q. 

In this work, we investigate the early growth of density fluctuations in spinodal region 
in charge asymmetric nuclear matter. For this purpose it is sufficient to consider the linear 
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response treatment of dynamical evolution 8j. The small amplitude fluctuations of the 
single-particle density matrix around an equilibrium state p°) are determined by the 
linearized TDHF equations. The linearized TDHF equations for fluctuations of neutron and 
proton density matrices, Sp^(t) = p^(t) — are given by, 

i^5p X a (t) = Kh X a {t)] + [5U*(t),p° a ]. (5) 

Since for infinite matter, the equilibrium state and the associated mean-field Hamiltonian h° a 
are homogenous, it is suitable to analyze these equations in the plane wave representations, 
d 

ih— < pl\5p a (t)\p 2 >= 

[ea(Pl) - £a0~2)] < Pl\8Pa{t)\P2 > [pa{fl) ~ Paifl)} < Pl\5U a (t)\f 2 > . (6) 

According to the basic assumption, matrix elements of the initial density matrix are Gaus- 
sian random numbers. In the plane wave representation the second moments of the initial 
correlations is given by, 



< pl|*Pa(0)|pa >< p 2 '|<W0)bl' > = (7) 

5 ab {27ch) 6 5{p~[ - Pl)S{p 2 - f 2 ')- [Pa(fl)(l - Pai.fl)) + Pa(P2)(l ~ Paifl))}, 

where the factor (2irh) 6 arises from normalization of the plane waves. 



III. GROWTH OF DENSITY FLUCTUATIONS 



A. Spinodal Instabilities 

In this section, we apply the stochastic mean-field approach in small amplitude limit to 
investigate spinodal instabilities in charge asymmetric nuclear matter [11]. We note that 
the following quantity 

5h a (k,t) = 2 / j-Z- < p + hk/2\S Pa (t)\p- hk/2 > (8) 

defines the Fourier transform of the local density fluctuations of neutrons and protons. In 
this expression and in other formulas in this section, we omit the event label A for clarity of 
notation. We can obtain solution of the linear response Eq. ([6]) by employing the standard 
method of one-sided Fourier transform in time, 

5h a (k,u)= dte iut Sn a (k,t). (9) 
Jo 
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After transformation, we obtain a set of coupled algebraic equations for the Fourier trans 



forms of fluctuating parts of local neutron and proton densities 

[1 + F™ Xn (k, uj)]5n n (k, u) + F^ X n(k, u)5h p (k, u) = tA n (k, uo) (10) 

and 

[1 + Fg P x P (k, uj)}5h p (k, u) + Fg n x P (k, uj)5h n (k, u) = iA p (k, u). (11) 

In these expressions, derivative of the mean-field potential U a (n n ,n p ) evaluated at the equi- 
librium density F^ = (dUb/dn a ) denotes the zero-order Landau parameters and Xa(k,u) 
is the Lindhard functions associated with neutron and proton distributions, 

ft \ o f°° d "P Pa(p-frk/2)- Pa (p + hk/2) 
Xa(k,u) = -2 — — — =- . (12) 

J-ooi 2 ^^) hu-p-hk/m 
The source terms A a (k,u) are determined by the initial conditions, 

A a (k, u) = 2h ^ <P+ h ^ 2 \ 6 PM[P- h k/2> ^ 

J-oo (2vr/?,) 3 huj — p-hk/m 

The solution of the coupled algebraic equations for Fourier transform of density fluctuations 
is given by, 

5fi ^ w) = , [1 + F pp Xp (k, oo)}A n (k, u) - F^xnjk, oo)A p (k, to) 

e(k, uj) 

and 

5fi ^ w) = , [1 + FrXnjk, u)}A p (k, u) - F pn X P (k, u)A n {k, u) 

s(k,u) 

where the quantity 

e(k, w) = 1 + F™Xn(k, w) + F™ Xp (k, u) + [F™F™ - F^F pn ]Xn(k, cu) Xp (k, u) (16) 

denotes the susceptibility. 

Time dependence of Fourier transform of density fluctuations 5h a (k,t) is determined by 
taking the inverse transformation of Eqs. (|T4|) and (fT5j) [3]. The inverse Fourier trans- 
formations in time can be calculated with the help of residue theorem, keeping only the 
growing and decaying collective poles we find, 

8n a (k, t) = 5n+(k)e +Tkt + 5n~(k)e- Tk \ (17) 
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where the initial amplitude of density fluctuations are given by 

ft^g) _ ( t 1 + jfffrpft <o)]Mk,u) ~ K P Xn(k, u;)A p (k, u>) (ig) 
1 de(k,u))/du) 



a;==F«r fe 

and 

[ [1 + K n Xn(K u^A^lu) - F^Xpjk, tMJM u) \ _ (ig) 

I de{k,u)/du L= Tl T k 

Growth and decay rates = ^fiT^ are determined from the dispersion relation e(k,u) = 0, 
i.e. from the roots of susceptibility. 

In numerical calculations we employ the same effective Skyrme potential as in reference 



U a (n n , n p ) = a(-)+B (-) ^ + C (-) r a + - DAn + D'An'r a (20) 

where n = n n + n p and n' = n n — n p are total and relative densities, and r a = +1 for neutrons 
and r a = —1 for protons. The parameters A = —356.8 MeV, B = +303.9 MeV, a = 1/6 
and D = +130.0 MeV fm 5 are adjusted to reproduce the saturation properties of symmetric 
nuclear matter: The binding energy 6q = 15.7 MeV/nucleon and zero pressure at the 
saturation density no = 0.16 fm~ 3 , compressibility modulus K = 201 MeV and the surface 
energy coefficient in the Weizsacker mass formula a sur f = 18.6 MeV [14]. Magnitude of the 



parameter D' = +34 MeV fm 5 is close to magnitude given by the SkM* interaction 15]. 
The potential symmetry energy coefficient is C(n) = Ci — C^in/no) 01 with C\ = +124.9 MeV 
and C2 = 93.5 MeV. These parameters for the symmetry energy coefficient in Weizsacker 
mass formula, at saturation density gives a sym = £p(no)/3 + C(n Q )/2 = 36.9/3 + 31.4/2 = 
28.0 MeV. 

As an example, Fig. HJa) shows the growth rates of unstable modes as a function of wave 
number in the spinodal region corresponding to initial density n = 0.2 n and n = 0.4 n 
for initial asymmetry / = 0.0 at a temperature T = 5 MeV. The initial charge asymmetry 
is defined according to I — (n° n — n^)/{n^ + n p ). In this figure and also in other figures, 
solid-lines and dashed-lines show quantal and semi-classical results, respectively. Since, at 
low densities, wave numbers of most unstable modes are comparable to Fermi momentum, 
long-wavelength expansion of the Linhard function is not valid, and hence there is important 
quantal effect in the dispersion relation. At the initial density n = 0.2 n and the initial 
asymmetry / = 0.0, in the quantal calculations unstable modes are confined to a narrower 
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FIG. 1: Growth rates of unstable modes as a function of wave number in spinodal region corre- 
sponding initial densities and at a temperature T = 5 MeV. (a) for initial asymmetry / = 0.0 , 
(b) for initial asymmetry / = 0.5. 

range centered around wavelengths A ~ 8 — 10 fm, as compared to a broader range centered 
around A ~ 7fm in the semi-classical calculations. Growth rates in semi-classical framework 
are determined by the roots of semi-classical susceptibility, which is defined as in Eq. (|T6l) 
by taking the Lindhard functions Xa(k,ui) in the long wavelength limit given by Eq. ( |28l ). 
As a result, in the quantal calculations, the source has a tendency to break up into larger 
fragments as compared to the semi-classical calculations. Also, due to quantum effects, the 
maximum of dispersion relation is reduced by about a factor 3/4. Therefore, fluctuations 
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take more time to develop when quantum effects are introduced. At higher initial density 
n = 0.4 no , in both quantal and semi-classical calculations, dispersion relation is shifted 
towards longer wavelengths and it exhibits a similar trend as the one at the initial density 
n = 0.2 Uq. This quantal effect in dispersion relation of unstable modes was pointed out 
in the case of symmetric matter in a previous publication [if]]. Charge asymmetric nuclear 
matter exhibits a similar behavior as seen from figure QJb), which shows dispersion relation 
corresponding to initial densities n = 0.2 n$ and n = 0.4 no for initial charge asymmetry 
/ = 0.5 at a temperature T = 5 MeV. Figs. [2](a) and[2](b) shows the boundary of spinodal 
region in density-temperature plane corresponding to initial charge asymmetries I = 0.0 and 
I = 0.5 at a temperature T = 5 MeV for the unstable modes with wavelengths A = 9 fm 
and A = 12 fm, respectively. It is seen that with increasing charge asymmetry, spinodal 
region shrinks to smaller size in both quantal and semi-classical calculations. Furthermore, 
unstable modes are quite suppressed by quantal effects as compared to the semi-classical 
results in both symmetric and asymmetric matter. Results of semi-classical calculation are 



in agreement with the result obtained in reference 



B. Growth of Density fluctuations 



In this section, we calculate early growth of local density fluctuations in charge asym- 
metric nuclear matter. Spectral intensity of density correlation function a a b(k,t) is related 
to the second moment of Fourier transform of density fluctuations according to, 



d ab (k,t)(2Tr) 3 5(k — k') = 6n a (k,t)6n b (-k f ,t). (21) 

We calculate the spectral functions using the solution (|T7I) and employing expression ([7j) for 
the initial correlations to find, 

dab (k, t ) = E+(k,iF k ) {e2Fkt + e „ 2Ffct) + 2E- b (k,iT k ) (22) 

| [ds(k, u;)/<9u;] w=i r fe | 2 | [de(k, u) / 'du\ u=iTk \ 2 

where quantities E^ b (k,iTk), a,b = n,p, are given by, 

EUl ^r k ) = 4ft 2 (l + F^ Xp fl^ + m\F^ Xn f /J , (23) 
E*(k, Xk) = 4ft 2 (l + F™ Xn fi; + Ah\Fr Xp ) 2 11 (24) 
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FIG. 2: Boundary of spinodal region in density-temperature plane corresponding to initial charge 
asymmetries I = 0.0 and I = 0.5 for the unstable mode: (a) with wavelength A = 9 fm , (b) with 
wavelength A = 12 fm. 

and 

EU^ ^k) = -4^(1 + F™ Xp )F pn X P 11 - 4ft 2 (l + F™ Xn )F^ X n 1* (25) 

with 

y (27rft) d (^,r fc ) 2 + (p • hk m) 2 } 2 



10 



Semi-classical limit of these expressions are obtained by replacing the integrals 1^ and 
Xa{k,w) with following expressions in the long wave-length limit, 

= / (o^m ^ — ,~ fc r/ ^ ^W 1 ~ P*w (27) 
7 (27rfr) d (M\) 2 + (p • hk/m) 2 } 2 



and 



. rf 3 p (g- fifc/m) 2 d 

- 2 / , (2^)3 (^ + ^ . ^ /m)2 ^- ( 28 ) 



Figs. E{a) and [3(b) shows spectral intensity a nn {k,t) of neutron-neutron density corre- 
lation function as function of wave number at times t = and t = 50 fm/c for density 
n = 0.4 no and the initial charge asymmetry / = 0.5 at temperature T = 1 MeV and 
T = 5 MeV^, respectively. 

As mentioned above, in all figures solid-lines and dashed-lines indicate quantal and semi- 
classical results, respectively. As seen, in particular at towards the high end of the wave 
number spectrum, considerable quantal effects are present at initial fluctuations. Quantum 
statistical effects in the initial fluctuations become even larger at smaller temperatures. In 
fact at zero temperature, since the quantities I^{sc) becomes zero, spectral functions vanish 
cfab{k,t) = 0. However, in quantal calculations spectral functions remains finite even at 
zero temperature, reflecting quantum zero point fluctuations of the local density. Looking 
at the results at t = 50 fm/c , we observe that largest growth occurs over the range of 
wave numbers corresponding to the range of dominant unstable modes. At T = 5 MeV, 
magnitude of fluctuations is about the same in both quantal and semi-classical calculations. 
At the lower temperature T = 1 MeV, magnitude of fluctuations in the most unstable 
range is nearly doubled in quantal calculations as compared to semi-classical calculations. 
Fig. [3j(c) shows spectral intensity a nn (k,t) as function of wave number at times t = 
and t = 50 fm/c at a lower density n = 0.2 no for initial charge asymmetry / = 0.5 and 
temperature T = 5 MeV. At the lower density, growth rates of dominant modes in the semi- 
classical limit are considerably larger than those of quantal calculations. Consequently, the 
result of semi-classical calculations at time t = 50 fm/c overshoots the result of quantal 
calculations over the range of dominant modes. Fig. H] illustrates that the spectral intensity 
for symmetric matter has similar properties as for asymmetric matter with / = 0.0. 

We note that quantal effects enter into the spectral density in two different ways: (i) 
quantal effects in growth rates of modes and (ii) quantum statistical effects on the initial 
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FIG. 3: Spectral intensity a nn (k,t) of neutron-neutron density correlation function as function of 
wave number k at times t = and t = 50 fm/c for the initial charge asymmetry / = 0.5 : (a) 
for density n = 0.4 no at temperature T = 1 MeV, (b) for density n = 0.4 no at temperature 
T = 5 MeV , (c) for density n = 0.2 no at temperature T = 5 MeV . 

density fluctuations, which becomes increasingly more important at lower temperatures. We 
also note that in determining time evolution of 5n(k, t) with the help of residue theorem, 
there are other contributions arising from non- collective poles of susceptibility e(k, u) and 
from poles of A a (k,u). These contributions, in particular towards short wavelengths, are 
important at the initial state, however they damp out in a short time interval Therefore 
the approximate expression fl22|) for the spectral intensity a(k,t) of density fluctuations 
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FIG. 4: Same as Fig. 3 but for asymmetry / = 0.0. 

becomes more accurate for increasing time. 

Local density fluctuations 5n a (f,t) are determined by the Fourier transform of 5h a (k,t). 
In terms of spectral intensity a a b(k,t), which is defined in Eq. ( )2T|) . equal time density 
correlation function as a function of distance between two space locations is expressed as, 

/d 3 k 
e ifc.(f-f )a ab (k,t). (29) 

Total density correlation function is given by sum over neutrons and protons and cross- 
term, a(\f — r'\,t) = cr n n(\r — r'\,t) + u pp (\r — r '\,t) + 2a np (\f — r '\,t). The behavior of 
density correlation function as a function of initial density and temperature carries valuable 
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information about the unstable dynamics of the matter in the spinodal region. As an 
example, Figs. 0(a) and 0(b) illustrate total density correlation function as a function of 
distance between two space points at times t = and t = 50 fm/c at density n = 0.4 no 
and the initial charge asymmetry / = 0.5 for temperatures T = 1 MeV and T = 5 MeV, 
respectively. At temperature T = 5 MeV, quantal effects are not important, and hence semi- 
classical calculations provide good approximation for density correlation function. However, 
at lower temperature T = 1 MeV, semi-classical calculations severely underestimates peak 
value of density correlation function. Fig. [5]^c) shows density correlation function at times 
t = and t = 50 fm/c at a lower density n = 0.2 no for initial charge asymmetry 
/ = 0.5 and a temperature T = 5 MeV . On the other hand, at lower density, semi-classical 
approximation overestimates the peak value of the correlation function. As indicated above, 
this is due to the fact that growth rates of dominant modes in semi-classical limit are 
considerably larger than those obtained in quantal calculations. For asymmetry / = 0.0, as 
seen from Fig. El behavior of density correlation function is similar to the charge asymmetric 
case. Complementary to the dispersion relation, correlation length of density fluctuations 
provides an additional measure for the average size of primary fragmentation pattern. We 
can estimate the correlation length of density fluctuations as the width of correlation function 
at half maximum. Correlation length depends on density, and to some extend, depends on 
temperature as well. From these figures, we can estimate that the correlation length of 
density fluctuations is about 3.5 fin at density n = 0.4 no, and about 3.0 fm at density 
n = 0.2 no- 

During spinodal decomposition, initial charge asymmetry shifts towards symmetry in 
liquid phase while gas phase moves toward further asymmetry. As a result, produced frag- 
ments are more symmetric than the charge asymmetry of the source. This interesting fact 
is experimentally observed and it may provide a useful guidance to gain information about 
symmetry energy in low density nuclear matter. For each event, we can define perturbation 
charge asymmetry during early evolution of density fluctuations as, 



We are interested in the ensemble average value of this quantity, which can approximately 
be evaluated according to 



P t - 



5n n (f,t) - 5n p (f,t) _ [5n n {f,t)f - [Sn p (f,t)} 2 
Sn n (f, t) + 5n p (f, t) [Sn n (f, t) + 5n p (f, t)] 2 



(30) 



I 



(31) 



pt 



(t) + 2a np (t) + a pp (t)' 
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FIG. 5: Density correlation function a(x,t) as a function of distance x =| f — r ' | between two 
space points at times t = and t = 50 fm/c and the initial charge asymmetry I = 0.5 : (a) 
for density n = 0.4 no at temperature T = 1 MeV, (b) for density n = 0.4 uq at temperature 
T = 5 MeV, (c) for density n = 0.2 no at temperature T = 5 MeV. 

where cr a fe(t) = cr a6 (| r — r ' |= 0,t). The average value of the perturbation asymmetry is 
nearly independent of time. As an example, Fig. [7] shows this quantity as function of initial 
asymmetry at temperature T = 5 MeV for densities n = 0.2 n and n = 0.4 n . As a result 
of the driving force of symmetry energy, perturbation asymmetry drifts towards symmetry. 
At this temperature quantal effects do not play an important role and these calculations are 
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FIG. 6: Same as Fig. 5 but for asymmetry I = 0.0. 



consistent with results of ref. 11]. 



IV. CONCLUSIONS 



Recently proposed stochastic mean-field theory incorporates both one-body dissipation 
and fluctuation mechanisms in a manner consistent with quantal fluctuation-dissipation the- 
orem of non-equilibrium statistical mechanics. Therefore, this approach provides a powerful 
tool for microscopic description of low energy nuclear processes in which two-body dissi- 
pation and fluctuation mechanisms do not play important role. The low energy processes 
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FIG. 7: Perturbation asymmetry as function of initial asymmetry at temperature T = 5 MeV for 
densities n = 0.2 no and n = 0.4 no. 

include induced fission, heavy-ion fusion near barrier energies, spinodal decomposition of 
nuclear matter and nuclear multi-fragmentations. In this work we investigate quantal ef- 
fects on spinodal instabilities and early growth of density fluctuations in charge asymmetric 
nuclear matter. For this purpose it is sufficient to consider the linear response treatment of 
the stochastic mean-field approach. Retaining only growing and decaying collective modes, 
it is possible to calculate time evolution of spectral intensity of density correlation function 
and the density correlation function itself including quantum statistical effects. Growth 
rates of unstable collective modes are determined from a quantal dispersion relation, i.e. 
from the roots of susceptibility. Due to quantal effects, growth rates of unstable modes, in 
particular with wave numbers larger than the Fermi momentum, are strongly suppressed. 
As a result, dominant collective modes are shifted to longer wavelengths than those obtained 
in the semi-classical description with the same effective interaction, in both symmetric and 
asymmetric matter. The size of spinodal zone associated with these modes is reduced by 
the quantal effects. In calculation of density correlation function, quantal effects enter into 
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the description through the growth rates of the modes and through the initial density fluc- 
tuations. Quantum statistical influence on density correlation functions grows larger at 
lower temperatures and also at lower densities. Quantal effects appear to be important 
for a quantitative description of spinodal instabilities and growth of density fluctuations in 
an expanding nuclear system. Stochastic mean-field approach incorporates both one-body 
dissipation and fluctuations mechanisms in a manner consistent with dissipation-fluctuation 
theorem. Therefore, it will be very interesting to investigate spinodal decomposition of an 
expanding nuclear system in this framework. We also note that numerical effort in simu- 
lation of stochastic mean-field approach is not so much greater than the effort required in 
solving ordinary three dimensional time dependent Hartree-Fock equations. 
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